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We present the quantum phase diagram of the one-dimensional t- J z model for arbitrary spin (in- 
teger or half- integer) and sign of the spin-spin interaction J z , using an exact mapping to a spinless 
fcrmion model that can be solved exactly using the Bethe ansatz. We discuss its superconducting 
phase as a function of hole doping v. Motivated by the new paradigm of high temperature super- 
conductivity, the stripe phase, we also consider the effect the antiferromagnetic background has on 
the t-J z chain intended to mimic the stripe segments. 
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Introduction. Phase diagrams of generic models of 
strongly interacting quantum particles are considered 
fundamental to understanding the complex physical be- 
havior of cuprate superconductors, heavy fermion, and 
related compounds. It is rare to encounter situations 
where unambiguously these diagrams can be completely 
determined and only a few exceptional cases are exactly 
solvable. It is a purpose of this paper to show that the 
t-J z chain belongs to this latter class of models. 

A new paradigm in superconductivity springs up as 
a consequence of the growing body of experimental ev- 
idence suggesting that the quantum state of high tem- 
perature superconductors is a striped phase. Unlike con- 
ventional conductors where the charge carriers distribute 
in an spatially homogeneous way, the stripe paradigm 
assumes that carriers cluster into quasi one-dimensional 
(Id) channels. These channels act as domain walls sep- 
arating different antiferromagnetic (AF) domains. It is 
remarkable that experiments are consistent with a spin 
ordering that is 7r-shifted across the wall 0], indicating 
the topological character of these extended defects || . 

Motivated by this new paradigm Ref. || argued that 
planar models, with appropriate inhomogeneous mag- 
netic terms, breaking translational and local spin SU (2) 
symmetries are appropriate to understand neutron scat- 
tering and angle-resolved photoemission spectroscopy ex- 
periments in cuprates |Q| . It is interesting to understand 
why spin anisotropies are relevant to obtain substantial 
pair hole binding and whether the stripes themselves have 
important superconducting fluctuations. The simplest 
representation of a stripe segment is realized by a t-J z 
chain model. 

In this Letter we study the quantum phase diagram of 
the t-J z chain for arbitrary spin and sign of J z by us- 
ing an exact mapping to an attractive spinless fermion 
model, and solve this problem by using the Bethe ansatz 
integral equations. We then consider the effect of the AF 
boundaries on the stripe as an effective confining poten- 
tial and determine the resulting phase diagram. While 
a superconducting phase exists in both cases, the super- 
conducting region is more prominent in the latter. 

Model Hamiltonian. The Hamiltonian representing the 



Id t-J z model with L sites (equal to the length of the 
chain, i.e., lattice constant a—1) and M holes with open 
boundary conditions (BC) g] (the thermodynamic limit, 
L, M — ► oo with v = M/L finite, is performed at the end 
of the calculation), for arbitrary half-integer spin S, is 
H = f + H Jz with 

L-l 

T = — t T a ,a i T aj(T — c* a(T c a+1(T + H.c. , 

a=l 
<t£[-S,S] 



L-l 



H Jz - J z ^2 Sa S a+l > S a=/. a 
a=l o£[-S,S] 



Here, (^(c^) creates(annihilates) a fermion of z spin- 
component a in a Wannier orbital centered at a. The 
Hilbert space of states corresponds to a constrained space 
with no doubly occupied sites || . 

Consider the set of parent states, labeled by the string 
configuration cr , with M holes and L — M quantum par- 
ticles, | $o (<*)), defined as 

|$q(<t)) = \ o- 1 cr 2 Q'3 ■■ ■ o o o o o ■ ■ ■ ) , (1) 

L — M M 

where a a indicates the z-component of the spin of 
the particle at site a. These states are eigenstates of 
Hjz with energy E M (&) = J z E„=?' 1 o- a a a+1 , and z- 
component of the total spin S z = X^a^f 1 °a- 

^From a given parent state one can generate a sub- 
space of the Hilbert space M{ff) by applying the hopping 
operators T a a to the parent state and its descendants, 

l*i(*)>, 



!$!(<?)> =r L _ 



l*o(<?)> 



or, in general, 



=f a , u 



(2) 



(3) 



The dimension T> of the subspace SA(<f) is More- 
over, these different subspaces are orthogonal and are not 
mixed by the Hamiltonian H, although they can be de- 
generate. In the following we will only consider the AF 
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case J z > 0. At the end we will discuss two important 
generalizations: the ferromagnetic (FM) J z < and the 
arbitrary integer spin hard-core boson cases. 

Among all possible initial configurations the one cor- 
responding to the Neel string <r (i.e., <j a = (— 1) Q S), 
which is two-fold degenerate, turns out to be special. 
We want to show now that for a given number of holes 
M the subspace generated by the Neel parent state, 
A4(<xq) = Mo, contains the ground state. To this end, 
one has to realize that the kinetic energy matrix ele- 
ments ($i(<x)|T|<I>j(<x)) are the same for the different sub- 
spaces M.. Nonetheless, the magnetic matrix elements 
(§i(a)\Hj z \$j((T)) = 5ijAi(a) are different for the dif- 
ferent subspaces, with Ai(cr ) < Ai(cr) [Q. Notice that if 
one assigns o a = to the presence of a hole at site a, then 
Ai(a) = ^zEo=i cr a (7 a+i- Therefore, the Hamiltonian 
matrices Hfj (of dimension V x V) in each subspace M. , 
consists of identical off-diagonal matrix elements (H^ = 

Hfj , i 7^ i) and different diagonal ones. These hermi- 
tian matrices can be ordered according to the increasing 
value of the energy Em of their parent states (for fixed 
L and M). For any E M (a) < E M {a), H M ' — H M + B, 
where B is a positive semidefinite matrix. Then, the 
monotonicity theorem tells us that 



as XXZ model) with Hamiltonian [fu 



E k {a)<E k {a') Vfc = V 



(4) 



where Ek(&)'s are the eigenvalues of H M arranged in 
increasing order (E^er) < Ek+i{cr))- Therefore, we con- 
clude that the lowest eigenvalue of H must be in A4q, is 
Ei(cto), and is two-fold degenerate. 

Spinless Fermion Mapping. The next step consists in 
showing, within the ground state subspace A4q, that the 
Hamiltonian H maps into an attractive spinless fermion 
model. If one makes the following identification 



i in ■■■ oooo 

L—M M 



L—M 



M 



(5) 



i.e., any spin {c^ aS or c* a _ s ) maps into a single spinless 
fermion (&JJ in A^Oj it is straightforward to show that 
all matrix elements of H^ are identical to the matrix 
elements of the interacting quantum lattice gas 



L-l 



L-l 



H = -t (bib a+1 + H.c.) - J Z S 2 ^2 n a.n a+ i 



(G) 



in the corresponding new basis. In Eq. |^, n a = b\ c b a . 

Quantum Phase Diagram. The attractive spinless 
fermion model of Eq. ^| certainly has a superconduct- 
ing phase (i.e., correlation exponent K p > 1) ||. For 
arbitrary values of J Z S 2 , t, and hole density v, the spin- 
less model is equivalent (via a Jordan- Wigner transfor- 
mation) to a Heisenberg-Ising spin-i chain (also known 



E 

a=l 



J±Ksl +1 + sls y a+1 ) + J,, (s z a s z a+1 + s z a ) 



(7) 



and J_l = 2i, J|| — —J Z S 2 . In this new representation the 
spin up(down) density is v(l — v). Jj_ = represents the 
classical Ising limit while Jii = is the extreme quantum 
limit (XY model). In general, the exchange anisotropy 
parameter A = J||/./± < determines the physical na- 
ture of the correlations while J± defines the energy scale. 

There is a vast literature on this model that was ex- 
actly solved by the Bethe ansatz For |A| < 1, solu- 
tions of this model belong to the universality class called 
"Luttinger liquids" |l2| , with correlation functions char- 
acterized by power laws with non-universal exponents 
continuously depending on A. The correlation exponent 



K p is determined from 



K n = —v KV n 



v 2 nD r 



(8) 



where the Drude weight (or charge stiffness) D c is related 
to the velocity of charge excitations v p by D c = v p K p /tt 
with D c = i5 2 e($)/9($/i) 2 as usual |§, and k is the 
isothermal compressibility. On the other hand, k can be 
calculated from the variation of the ground state energy 
per site e{v) as {v 2 ^ 1 — d 2 ejdv 2 . 

At v = 1/2, several quantities and properties are 
known in closed analytic form. There is a Mott tran- 
sition at A = 1 (Umklapp scattering becomes relevant at 
A = 1, while it is irrelevant for |A| < 1). Moreover, the 
exact expressions for K p and v p can be determined from 
the Bethe ansatz [Ol 



irt sin /i 



(9) 



which implies D c = irt sin /i/[4/i(7r — (i)] with cos/i = A. 
The energy per site (|A| < 1) is given by 



, , s At Z" 00 , sinh(7r- u)x 

e(l/2) = — - 2tsin^ / dx — i — r^— 

2 Jo cosh fix smh nx 



(10) 



Thus, one obtains \ < K p < 1 in the region > A > 
— 1/ V2, and K p > 1 (superconducting correlations dom- 
inate at large distances) for —1 < A < — l/y/2. At 
A = — 1, there is a transition to a phase segregated state 
(k = 2^/[7rf sin \i (ir — /i)] diverges). 

For A = 0, the system reduces to a free spinless 
fermion system with energy per site e = — — sixifnu), 
stiffness D c = — e/2, and k^ 1 = —n 2 v 2 e. This trivially 
corresponds to K p — i . Also the cases v — > and v — > 1 
map to free fermions independently of the value of A, 
therefore, K p = i. The value A = — 1 is also special: 
After the unitary transformation s'^' v — > (— l) a s^ ,y , the 
Hamiltonian H xxz maps into the FM Heisenberg model in 
a magnetic field Jn (full SU(2) symmetry is recovered). 
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Here e = — 5 + (5 — v \ implying the opening of a gap 
with a diverging k, i.e., A — — 1 determines the line of 
phase segregation for all densities v. 



(a) 
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FIG. 1. Quantum phase diagrams of the (a) t- J z chain, and 
(b) modified t-J z model which includes the effects of the AF 
background where the chain is embedded. There are three 
different quantum phases: phase segregated (PS), supercon- 
ducting (SC), and metallic (C) phases. The last two be- 
long to the "Luttinger liquid" universality class. Points with 
p = 1 known in closed analytic form are: yi = (- 

(0, 2-V2), and y 3 = (cos fvr, 



K 



'73' 2; 
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Away from v = 1/2 and the special limiting cases 
discussed above, the quantities K p and v p are obtained 
from the numerical solution of the Bethe ansatz inte- 
gral equations (TJJl^]. To calculate v p one needs to de- 
termine hole excitations with a well-defined momentum 
q and energy Ae with respect to the ground state e. 
We find the velocity of this elementary excitation from 
v p = ]im q —>odAe/dq, and together with the numerical 
second order derivative of e(v), we determine the cor- 
relation exponent K p . For |A| < 1 the excitations are 
gapless. The resulting quantum phase diagram is shown 
in Fig. |l|(a). Notice that the largest superconducting 
region corresponds to v = 1/2. 

It is interesting to determine the influence of the an- 



tiphase domain wall (ADW), associated with each charge, 
on the spin-spin correlations of the metallic and the su- 
perconducting phases. It is well-known |L6| that the 
charge structure factor of the spinless model has a peak 
at k = 2kp = 2ttv. Since each charge of the effective 
model carries an ADW, the spin structure factor, S(fc), 
will peak at k — 7r ± 2-kv in the metallic phase. In the 
superconducting phase a broad peak at k — it is obtained 
for S(fe) since the pairs do not carry an ADW. 

Effect of AF Boundaries on the Chain. It is also inter- 
esting to study the effect of the AF background in which 
our stripe segments are embedded. This background pro- 
vides a strong BC that results in an additional attractive 
(confining) potential for the holes on the stripe. In this 
way, an enhanced superconducting region is expected. In 
fact, the influence of the AF background on the stripe is 
equivalent to the effect a staggered magnetic field (STM) 
B s (see Fig. ||). Since each hole carries an ADW, the 
staggered field gives rise to a confining linear potential 
between the a and the a + 1 holes for even a. There- 
fore, the effect of a STM is to bind pairs of holes tightly 
by a string potential. These pairs interact as hard core 
bosons. In the very dilute limit one expects these bosons 
to condense at T = for any value of B s /t and J z /t 
lower than some critical value which gives rise to phase 
segregation. 

In the limit B s /t ^> 1 the model can be solved analyt- 
ically for any concentration of hole pairs. In this limit 
the problem reduces to nearest-neighbors (NN) pairs 
of holes moving into an AF background (see Fig. 2). 
Each pair can hop to its NN with an effective hopping 
t c ff = 2t 2 /(B s S + J Z S 2 ). In addition, there is an at- 
tractive J Z S 2 interaction between NN pairs which comes 
from the second term of Eq. ^ If we map each pair into 
a spinless particle and each spin into an empty site |l7| 



|ooTIT|oot| 



•oooo»oo 



(L,M) 



(L,M) 



(11) 



the problem reduces to the spinless Hamiltonian Eq. ||. 
But as each pair is replaced by an effective particle, the 
number of particles and the length of the chain for the 
effective spinless problem are 



L = L 



M 

T 



M 



M 

Y 



(12) 



where v = M/L — v/{2 — v). Here, as in the FM case, 
we can use closed or open BC. The sign that arises af- 
ter a cyclic permutation of fermions is absorbed in the 
BC. For an odd(even) number of fermions these BC are 
periodic (anti-periodic). This corresponds to M = An + 2 
[M = An) of the original problem. The relations between 
energies and charge velocities are [{18[ 



i{y) = e{y) (l - ~) , v p {v) = v p {v)- 



(13) 
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where e(v) and v p {y) are the energy per site and charge 
velocity of the corresponding spinless model of concentra- 
tion v. Therefore, simple algebraic manipulations lead to 

K p (u) = K p {v) (2 - vf , (14) 

and the phase diagram is depicted in Fig. |l|(b). 



J J 




J z J z 

FIG. 2. Schematic representation of the effect of an AF 
background on a stripe segment. In this picture we assume 
that the concentration of holes in the stripe is such that 
the there is no 7r shift between AF domains. If there were 
anti-phase domains, then, there is no confining string poten- 
tial even though two holes still like to be in adjacent sites. 

For completeness, we would like to mention that the 
mapping of the low energy spectra of the t-J z model into 
the spinless Hamiltonian, Eq. |^, is also valid for the FM 
case, i.e., J z < 0. In this case, the magnetic background, 
which is replaced by empty sites in the spinless model, is 
FM. Notice, however, that the effective spinless model is 
also attractive (A < 0). This implies that the dynamics 
of the charge degrees of freedom in an AF background 
is the same as in the FM one. But in the latter case, 
the charges do not carry an ADW. Moreover, the map- 
ping does not depend upon the statistics of the quantum 
particles. In other words, we could also apply these con- 
cepts to constrained quantum particles with integer spin 
S 1 , i.e., hard-core bosons. In the large spin S limit, the 
quantum phase diagram of the t-J z model approaches 
the one of isotropic t-J Hamiltonian. Notice, however, 
the qualitative similarity between the phase diagram in 
Fig. ^](a) and the one for the isotropic spin-i t-J model 
obtained numerically |1S[] . Finally, the solution can be 
trivially extended to the t-J z -V model, where V repre- 
sents a NN density-density interaction. The effect of V is 
simply to renormalize the spinless fermion interaction in 
Eq. |^. Furthermore, it is simple to prove that there is a 
family of bilinear- biquadratic spin-1 chain Hamiltonians 
that can be mapped onto a t-J z -V model and, therefore, 
its low energy physics is exactly solvable [E0| . 

In summary, we presented the exact quantum phase di- 
agram of the t-J z chain model for arbitrary spin S, parti- 



cle statistics, and sign of the magnetic interaction J z . We 
also exactly determined the phase diagram of a modified 
t-J z chain that includes the effects of a strong antiferro- 
magnetic background. A metallic, superconducting and 
segregated phases characterize these two phase diagrams. 

We thank A. A. Aligia, S. Sondhi, and J.E. Gubernatis 
for useful discussions, and A. A. Aligia for helping us with 
the Bethe ansatz equations. Work at Los Alamos is spon- 
sored by the US DOE under contract W-7405-ENG-36. 
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